NEGATIVE REFLECTION IN A FARADAY CHIRAL MEDIUM 
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ABSTRACT: The four wavenumbers associated with planewave propagation in a Faraday 
chiral medium (FCM) with relatively huge magnetoelectric coupling give rise to enhanced pos- 
sibilities for negative-phase-velocity propagation and therefore negative refraction. They can 
also give rise to the phenomenon of negative reflection. In particular, for a nondissipative ex- 
ample, we deduced that an incident plane wave with positive/negative phase velocity can result 
in a negatively reflected plane wave with negative/positive phase velocity, as well as a positively 
reflected plane wave with positive/negative phase velocity. 
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1. INTRODUCTION 

The scope for exotic — and potentially useful — electromagnetic phenomenons to arise is 
greatly extended in anisotropic mediums, and especially bianisotropic mediums, as compared 
with isotropic mediums [1, 2]. This letter concerns a particular bianisotropic medium called a 
Faraday chiral medium (FCM). FCMs combine natural optical activity (as exhibited by isotropic 
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chiral mediums [3]) with Faraday rotation (as exhibited by gyrotropic mediums [4, 5]). They may 
be theoretically conceptualized as homogenized composite mediums arising from the blending 

together of isotropic chiral mediums with either magnetically biased ferrites [6] or magnetically 
biased plasmas [7]. Accordingly, a wide range of constitutive parameters for FCMs may be 
envisaged. 

In a recent study, we identified the propensity of FCMs to support planewave propagation 
with negative phase velocity, particularly when the magnetoelectric constitutive parameters 
are relative large [8]. Similar behavior in isotropic chiral mediums and their nonreciprocal 
counterparts can be traced back to 1986 [9], but is only nowadays being considered carefully 
[10, 11], following the successful exhibition of the phenomenon of negative refraction by isotropic 
achiral materials [12]. In the following sections we consider the phenomenon of negative reflection 
in FCMs. Negative reflection has previously been described for isotropic chiral mediums [13, 14] 
and anisotropic dielectric mediums [15], but not to our knowledge for bianisotropic mediums. 

In the following sections, an exp(— iwt) time dependence is assumed; the permittivity and 
permeability of free space are denoted by eo and /xq, respectively; 3-vectors are underlined, with 
the caret (") denoting a unit vector; and 3x3 dyadics are double underlined. 

2. PLANEWAVE PROPAGATION 

In preparation for solving a boundary-value problem in Section 3, let us consider planewave 
propagation in a FCM characterized by the frequency-domain constitutive relations [16, 17] 

D{r) = i • E{r) + C * H{r) 
Biz) = • E{r) + jj, • H{r) 

with constitutive dyadics 
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Simplifying for the sole purpose of exemplifying negative reflection, we choose the distinguished 



axis of the FCM to coincide with the z axis; i.e., u = z. 

Let us also confine ourselves to propagation in the xz plane. Plane waves with field phasors 

E_{r) = AE_Q exp (ik • r) 

BAr) = AI^Q exp (ik • r] 
can propagate in the FCM, where A is a complex-valued amplitude and the wave vector 
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k = KX + ^/l^^—K^z_ = kk = k{xsm^/; + icos^p) (4) 

contains the real-valued quantity k that is fixed by the incidence conditions in Section 3. The 
angle is defined in the clockwise sense with respect to the +z axis in the xz plane. 

In order to determine k, E^q, and ^q, we substitute (3) in the Maxwell curl postulates 
V X -E(r) = iujB{r) and V x H_{r) = —iujD{r)^ and then use the constitutive relations of the 
FCM. In general, we find four independent wavevectors for a specific n — namely. 
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hj = + y /c| — z = kj{xsmtjjj + zcos i/jj) , j G [1, 4] , (5) 

Corresponding to each kj, there exist E^qj and I^qj = • (l/uj) {kj x /) + ^ • E^qj. Thus, 
for the purposes of Section 3, the complete representation of plane waves in the FCM for a 
specified k is 

^(i:) = Ei=i^i^o, exp (tk^-r) 

wherein we have neglected the remote possibility of the existence of Voigt waves [18]. 
The associated time-averaged Poynting vector is given by 

4 4 

P{r) = E Pjir) = 2 E I^jI' i^oj X ^o, ) exp [-2Im {k^) • r] . (7) 

The phase velocity of the j*'* plane wave in (6) is positive if Re (^kj) • Pj{r) > and negative if 
^e{kj) -Pjir) < 0. 

3. BOUNDARY-VALUE PROBLEM 

Suppose that the half-space z > is filled with a FCM described by the constitutive relations 
(1), whereas the half-space z < is vacuous. A distant source in the FCM-filled half-space 



supposedly launches a plane wave towards the boundary z = 0. The wave vector of this plane 
wave lies wholly in the xz plane, and its projection on the boundary z = is denoted by k. 

Equation (6) holds in the region z >0 with the following stipulations: 

(i) The indexes j = I and j = 2 are reserved for plane waves that transport energy towards 
the boundary z = 0, i.e., 



These two indexes thus identify incident plane waves. 

(ii) The indexes j = 3 and j = 4 are reserved for plane waves that transport energy away 
from the boundary z = 0, i.e., 



These two indexes thus identify reflected plane waves. 

(iii) Either Ai = or A2 = 0, but not both. When Ai 7^ 0, tpi can be called as the angle of 
incidence; when A2 ^ 0, is the angle of incidence. 

If wavevectors of the incident and a reflected plane wave are oriented in either (i) the same 
quadrant of the xz plane or (ii) diagonally opposite quadrants of the xz plane, then the reflection 
is called negative; otherwise, the reflection is called positive. 

In the vacuous region, the transmitted field phasors are given by 



z'Pj{r)<0, je[l,2]. 
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z'Pj{r)>0, iG[3,4]. 
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z<0, 



(10) 



H{l) = Vo (.BsP - Bpy) exp [i{KX - az)] 



where Bg and Bp are unknown amplitudes for the s and p polarization states, and 



p = {ax + Kz_) I ko 



a = 



2 



> 



(11) 
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Imposition of the usual boundary conditions across the boundary z = yields four algebraic 
equations: 

4 

x'J2^j^j = Bp(^/ko, (12) 

4 

y'^A,R,^ = Bs, (13) 

4 

x-^Aj Hoj = Bs a/ (kom) , (14) 

4 

y-^AjEo^ = -Bp/r,o. (15) 

We recall that either = or ^2 = 0, but not both, for our present purpose. 

4. A NUMERICAL ILLUSTRATION 

In order to provide a numerical illustration, we considered the nondissipative regime wherein 
p, pg, Pz & ^ {p = e,^,/x), and kj G M G [1,4]. For the results presented in this section, we 
fixed e = 3, eg = 1, = 2.5; ^ = 10, = 6, = 7; /x = 2, pg = 0.5 and pz = 1-5. 

The values of ijj for each value of k are plotted in Figure 1. Let us focus attention on a 
specific value of k; e.g., k = Ako. As highlighted in Figure 1, the values of ■0 at k = Ako are 
ipi = 83.25° (on thin solid curve), ip2 = 167.63° (on thin dashed curve), V'a = 39.78° (on thin 
dashed curve), and V'4 = 160.09° (on thin solid curve). 

The orientations of the wavevectors kj and the time-averaged Poynting vectors P_j at k = 4ko 
in the FCM are presented in Figure 2. This figure provides confirmation that the wavevectors 
with orientation angles '0i^2 are associated with energy flow towards the boundary z = 0, 
whereas the other two wavevectors are associated with energy flow away from that boundary. 
In addition, the plane waves indexed 1 and 4 have negative phase velocity, while the plane waves 
indexed 2 and 3 have positive phase velocity. 

Now, let us consider further the reflection of the incident plane wave indexed 1 by letting 
Ai = 1 and A2 = 0. This case results in a positively reflected plane wave with wavevector 
orientation angle 1/^4 and a negatively reflected plane wave with wavevector orientation angle 
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tps- The positively reflected plane wave has negative phase velocity whereas the negatively 
reflected plane wave has positive phase velocity. 

Next, we let Ai = and = 1 in order for the plane wave indexed 2 to be incident on 
the boundary z = 0. In this case, the reflected plane wave with wavevector orientation angle 
-03 is positively reflected and has positive phase velocity, whereas the plane wave indexed 4 is 
negatively reflected and has negative phase velocity. 

In order to check that the reflected plane waves described in Figure 2 have nonzero ampli- 
tude, the system of algebraic equations (12)-(15) was solved at k = Ako. With E_qj • E_*qj = 1 
(j G [1,4]), we obtained: = 0.006 - 0.079i, A4 = -0.892 - 0.372i, Bs = 0.078 - 0.306i, and 
Bp = -0.253 - 0.036i when ^1 = 1 and A2 = 0; and A3 = -0.921 - 0.477i, A4 = 0.073 - 0.071i, 
Bs = 0.106-0.380Z, and Bp = 0.329+0.070i when ^1 = and A2 = 1. Therefore, the negatively 
reflected waves do indeed have nonzero (but small) amplitude. We note that the refracted plane 
waves (in the vacuous region z < 0) are evanescent, since a = iko^/T5 at «; = Akg. 

The amplitudes of the negatively reflected plane waves increase if the vacuous region z <0 
is fllled by a perfect electric conductor. Then, we must set Bg = Bp = and the algebraic 
equations (12) and (13) have to be solved to find A3 and ^4. At k = 4fco with Eqj • E*Qj = 1 
{j G [1,4]), we found: A3 = -0.884 and A4 = -0.089 when Ai = 1 and A2 = 0; and A3 = 0.102 
and A4 = -1.139 when = and A2 = 1. 
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5. CONCLUDING REMARKS 

The four independent wavenumbers associated with planewave propagation in a FCM can 
give rise to a host of complex electromagnetic behavior. With a nondissipative example, we 
have shown that an incident plane wave with positive phase velocity can result in a negatively 
reflected plane wave with negative phase velocity, as well as a positively reflected plane wave 
with positive phase velocity. Also, an incident plane wave with negative phase velocity can result 
in a negatively reflected plane wave with positive phase velocity, as well as a positively reflected 
plane wave with negative phase velocity. Thus, negative reflection is characterized by a reversal 
in the co/contra directionality of the wavevector and the time-averaged Poynting vector. These 
negative reflection characteristics add to an already rich palette of electromagnetic responses 
supported by linear bianisotropic mediums [2]. 
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Figure 1: The wavevector orientation angle V (in degree) plotted against k (normalized with 
respect to ko). There are four possible values of "0 at each value of k. The vertical line indicates 
K = 4ko. 
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